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t~| , Abstract 

The phase structure of a four- and eight-fermion interaction model 
^\J , is investigated at finite temperature and chemical potential in arbitrary 

^ ' space-time dimensions, 2 < D < i. The effective potential and the gap 

r^ , equation are calculated in the leading order of the 1/N expansion. If 

the first order phase transition takes place, the phase boundary dividing 
the symmetric and the broken phase is modified by the eight-fermion 
interaction. 



o. 

S ■ 1 Introduction 

Models of particle physics can be classified based on its symmetry. Some 
/\ • apparent symmetry in the models is broken spontaneously. The basic in- 

^ . gredient of the standard electroweak theory is the spontaneous symmetry 

breaking of the SU{2) ® U{1) gauge symmetry and grand unified theories 
are constructed on the basis of the spontaneous breaking of higher gauge 
symmetry at high energy scale. The phase transition has played a decisive 
role in the scenario of cosmological and astrophysical evolution. There is a 
possibility to inspect models of particle physics in cosmological and astro- 
physical objects. It is quite interesting to make investigations on the spon- 
taneous symmetry breaking under the circumstance of the early universe, 
high temperature, chemical potential and so on. 

As is known, a strong gauge interaction between fermions is subject to 
chiral symmetry breaking in a dynamical way at vanishing temperature and 



chemical potential. The broken symmetry is restored at high temperature 
and chemical potential. The symmetry breaking takes place in the strong 
coupling region in the gauge field theories. Non-perturbative treatment is 
necessary to study its phase structure. We often use phenomenological ef- 
fective models to evaluate non-perturbative effects for the gauge theories. A 
four-fermion interaction model is one of the simplest models where the chiral 
symmetry is broken dynamically [U E] . It is also useful to describe meson 
properties as a low energy effective theory of QCD. 

In two space-time dimensions four-fermion interactions are given by marginal 
operators. Hence the model is renormalizable. In three space-time dimen- 
sions the model is solvable and renormalizable in the sense of the 1/A^ ex- 
pansion. We can also regard the model in the space-time dimensions less 
than four as a regularized model in four dimensions [3l H] . Many works have 
been done to study the phase structure of a four-fermion interaction model 
at finite temperature and chemical potential in a space-time dimension less 
than four. The critical temperature in the Gross-Neveu model has been 
found at the large A^ limit in Ref. [5j. The influence of chemical potential is 
discussed in Ref. [B]. The phase structure at finite temperature and chem- 
ical potential has been studied in two [Zl El [9], three [10], four [11] and an 
arbitrary space-time dimensions [H [T^ [13] . 

It should be noticed that there is no reason to omit higher dimensional 
operators in a phenomenological effective model at high temperature and 
chemical potential. For example, 't Hooft consider a six fermion interac- 
tion in a low energy effective model in QCD to deal with the axial anomaly 
|14j . An eight-fermion interaction has been discussed to describe the mass 
spectrum for the nonet pseudoscalar meson [151 IT6] . Higher derivative in- 
teraction has been studied in Ref. |17] . These interactions are given by 
irrelevant operators. The renormalized coupling constant should be diver- 
gent at a high energy scale. Thus a larger contribution is expected from 
irrelevant operators at high temperature and chemical potential [18l [19] . 

In the present paper we consider scalar type four- and eight-fermion in- 
teractions and investigate the contribution of the eight-fermion interaction 
to the phase structure at finite temperature and chemical potential in arbi- 
trary space-time dimensions, 2 < D < 4. In Sec. 2 we define a model with 
four- and eight-fermion interactions considered here. The effective potential 
is calculated in the leading order of the 1/A expansion. In Sec. 3 we intro- 
duce the temperature and chemical potential in the model according with 
the imaginary time formalism. In Sec. 4 we numerically evaluate the behav- 
ior of the effective potential with as the temperature and cheinical potential 
vary. Observing the minimum of the effective potential, we derive the phase 
boundary of the discrete chiral symmetry. In Sec. 5 we give some concluding 
remarks, 



2 Four- and Eight-Fermion Interaction Model 



Here we consider a model with only scalar type four- and eight-fermion 
interactions to find a characteristic feature of the contribution from a higher 
dimensional operator to dynamical symmetry breaking at finite temperature 
and chemical potential. We extend the Gross-Neveu type model to have 
a scalar type eight-fermion interaction in arbitrary space-time dimensions, 
2 < D < 4 and start with the action 
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where the lower index i represents the flavors of the fermion field %[), N is the 
number of the fermion species and Gi, G2 are coupling constants. Below we 
neglect the flavor index for simplicity. The action ([T]) is invariant under the 
discrete chiral transformation, 
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(2) 



This chiral symmetry prohibits the action to have a mass term. If the 
composite operator, ■0(j;)'0(j;), develops a non-vanishing expectation value, 
above chiral symmetry is eventually broken. 

In practical calculations it is more convenient to introduce auxiliary 
fields, s and a, and rewrite the action ([1]) with the equivalent one [2_1 J , 
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Replacing the auxiliary fields in the action ([3]) by the solution of the equa- 
tions of motion, 
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we reproduce the original action ([T]). The expectation value for the composite 
operator, ipip is measured by the one for the auxiliary field a. Thus we regard 
the expectation value, {a), as an order parameter for the chiral symmetry 
breaking. If the non- vanishing expectation value is assigned to s, there 



appears a mass term for the fermion field, ip. The dynamically generated 
fermion mass, md = (s), is given by a function of {a). 

We want to calculate the order parameter, (a), under the ground state. It 
is found by observing the minimum of the effective potential. In the leading 
order of the 1/A^ expansion the effective potential is written as [2^ 
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where S{x, x; m) is the spinor two-point function and tr denotes trace with 
respect to spinor index. This expression is also applicable to the state in 
thermal equilibrium. For vanishing temperature and chemical potential the 
two-point function is given by 

S{-,V:s)=l^e-^'^^-y)—^-, (7) 
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Substituting Eq.([7|) to Eq.dH]) and performing the integration, we obtain the 
effective potential for vanishing temperature and chemical potential. 
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where we put the suffix for Vo(s, a) to keep the memory that T = fj, = 0. 

On the ground state the expectation value of the auxiliary fields, s and 
a, should satisfy the stationary conditions. 
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Thus the auxiliary field, s, is described by a function of a, 
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(9) 
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(11) 



Using Eqs.(l9]) and (fTT]l . we obtain the self-consistency equation for the dy- 
namically generated fermion mass. 
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This equation has a trivial solution s = 0. Since T (1 — D/2) is negative for 
2 < D < 4, Eq ()12p has a real and non-trivial solution only if either Gi or 
G2 is negative. In the four-fermion interaction model, G2 = 0, the discrete 
chiral symmetry is broken for a negative Gi [12]. The solution is given by 
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On the stationary condition Eq. ()10p . the effective potential 
as a function of the auxiliary field, a. 



is written 
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where we normalize all the mass scale by uiq and set 
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We can regard the expectation value, {a)^ as the order parameter of the 
chiral symmetry breaking. On the ground state it is given by the minimum of 
the effective potential ([HI). The extremum of the effective potential satisfies 
the stationary condition, 

1 dVo{a) 
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This equation has a trivial solution. 



and non-trivial ones, 
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and 
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We plot a typical behavior of the non-trivial solutions ()18p and ()19p as a 
function of ^f in Fig. [TJ As is shown, the solution (|18p appears for a negative 
5 and disappears at the limit, G2 — ?• 0. In the present paper we assume 
that the discrete chiral symmetry is broken for vanishing temperature and 
chemical potential and take a negative value for Gi. Hence the eight-fermion 
coupling, G2, is positive definite for a negative g. The gap equation of the 
four-fermion interaction model is reproduced at the 5 — )• limit of Eq. ()19p . 

It should be noted that the stationary condition is only the necessary 
condition for the minimum of the effective potential. One of the solutions 
(fT7|) . (fT8]l and (fT9]l gives the expectation value for a on the ground state. 
The behavior of the effective potential (|14p is shown in Fig. [2] as ^f varies 
at D = 3. We clearly observe a local extremum which corresponds to the 
solution P^ . It disappears for g = —0.1 

Evaluating the effective potential, we find the the expectation value, {a), 
on the ground state. From Eq. (jlip the dynamical fermion mass, m^, is 
described as a function of the expectation value, (a), 
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In Fig. [3] we plot the dynamical fermion mass at -D = 3. We observe a 
small mass gap where the global minimum moves between the solutions (jlSp 
and (|19p . On the right-hand side of the gap the minimum of the effective 
potential is derived by Eq. (fT9]l . For the solution (fT8]l we find an simple 
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Figure 1: Solutions (fT8|) and (fT9]) as a function of g at D = 3. 




Figure 2: Behavior of the effective potential for g 
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Figure 3: Dynamical fermion mass as a function of g at D = 3. 



analytical expression for the dynamical fermion mass, 
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Thus the dynamical mass is given by Eq. (|21|) on the left-hand side of the 
mass gap in Fig. [3l 



3 Effective Potential at Finite Temperature and 
Chemical Potential 

To discuss the phase structure of the model we calculate the effective po- 
tential at finite temperature and chemical potential. According to the imag- 
inary time formalism, the spinor two-point function at finite temperature 
and chemical potential is obtained from Eq. ([7]) by the Wick rotation and 
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where /3 = l/{kBT) with /c^ the Boltzmann constant and T the tempera- 
ture, and fj, is the chemical potential. Thus the spinor two-point function is 
modified to be 

S{x,y;z) = i^y f ^-Z^e-'^^--y^—^ . (25) 

/5^y (27r)^-i Ji'kf.-s ^ ' 

Inserting Eq. (j25|) into Eq.Q, we obtain the effective potential at finite tem- 
perature and chemical potential, 
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We can perform the summation and integrate over angle variables in this 
equation. Then we obtain 

V{s,a) = Vo{s,a) + V^^{s), (27) 
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It should be noted that the thermal correction, V^'^{s), is positive-definite. 
Thus the temperature and chemical potential effect suppresses the chiral 
symmetry breaking. The numerical uncertainty increases at the zero-temperature 



limit of Eq. ()28p . To study the case at T = we take the hmit analyticahy 
and use the fohowing expression, 
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Integrating over space components in Eq. ()26p . we obtain another expression 
for the effective potential [12], 
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This expression is equivalent to the one in Eq. (|28p . We use it for the analytic 
study of the critical point. 

Since V^'^{s) is independent on the auxiliary field a, the stationary con- 
dition (jlOp is not modified by the thermal correction. Hence the effective 
potential ([27j) is described as a function of the auxiliary field, a, on the 
stationary condition (fTO|l . 






Via) = Voia) + V^^{s = a + g^]. (31) 



Below we evaluate the effective potential (|3ip at finite temperature and 
chemical potential and study the phase structure of the four- and eight- 
fermion interaction model. 

4 Phase Structure at Finite Temperature and Chem- 
ical Potential 

The broken chiral symmetry is restored at high temperature and chemical 
potential. In Ref. [12] the phase diagram is illustrated for the four-fermion 
interaction model in arbitrary dimensions, 2 < L) < 4. Here we extend the 
analysis to the four- and eight-fermion interaction model. 



4.1 Phase boundary for the second order phase transition 



Before the full analysis of the phase structure, we analytically consider a 
phase boundary in a restricted case. We assume that the phase transition 
is of the second order. If the broken chiral symmetry is restored through 
the second order phase transition, the order parameter, (a), continuously 
disappears at the critical point. The order parameter satisfies the stationary 
condition for the effective potential ()3ip . 
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For T = fi = the stationary condition has the solutions ()17p . (jlSp and ()19p . 
The constant solutions ()17p and (jlSp also satisfy the stationary condition 
p2p at finite temperature and chemical potential. The other solution ()19p 
is modified as 
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For the second order phase transition the order parameter develops contin- 
uously at the critical point. Thus the critical value for the temperature and 
chemical potential is obtained by taking the o" — t- limit in Eq. ()33p . 
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where ^ (z, a) is the generahzed zeta function. Taking the hmit fi 
Eq. ()34p we obtain the critical temperature for /i = 0, 
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It is clearly seen that Eqs. ([M|) . ([35]) and ([36]) are independent on the 
eight-fermion coupling constant. Therefore the critical values receive no cor- 
rection from the eight-fermion interaction for the second order phase tran- 
sition in the leading order of the \/N expansion. The contribution from 
the irrelevant operator vanishes at the critical point for the second order 
phase transition. Indeed, Eqs. ([M]) . ([35]) and ([36]) exactly reproduce the one 
obtained in the four-fermion interaction model 1121. 



4.2 Dynamical fermion mass 

In the broken phase the fermion mass is generated dynamically. In our 
model the fermion field, ifj, acquires a finite mass, if a non-vanishing value 
is assigned for {a). In the symmetric phase the fermion mass term is pro- 
hibited. Thus the phase boundary is found by observing the dynamically 
generated fermion mass. Here we numerically evaluate the effective poten- 
tial ([3ip and calculate the dynamically generated fermion mass as a function 
of the temperature and chemical potential. 

We want to understand the contribution from the eight-fermion interac- 
tion. For this purpose it may be instructive to review the behavior of the 
dynamically generated fermion mass in the four-fermion interaction model, 
5 = 0. As is shown in Ref. |12] . the fermion mass vanishes above the critical 
temperature or chemical potential. The first and the second order phase 
transition coexist for 2 < Z? < 3. Only the second order phase transition 
takes place for 3 < Z) < 4. 

In Figs, m [5] and [6] we plot the behavior of the dynamical fermion mass 
for 5 = 0. In the case of the first order phase transition a mass gap appears 
at the critical point. Hence we can clearly distinguish the order of the phase 
transition. For D = 2.5 we observe the mass gap for kT = 0, 0.2?7io and 
/i = 0.8?TT,o in Fig. Hj (a) and (b), respectively. The mass gap disappears 
for higher temperature. There is the first order phase transition at low 
temperature. In three dimensions the first order phase transition is observed 
only for T = 0. There is no first order phase transition for D = 3.5. In 



11 




0.4 0.8 1.2 1.6 

|i/mo 



1.4 

1.2 

1 

I 0.8 

1 0.6 

0.4 

0.2 






)i=0.4mi, 



0.2 0.4 0.6 0.8 1 1.2 

kT/mo 



(a) kT /rriQ is fixed. 



(b) fi/rriQ is fixed. 



Figure 4: Dynamical fermion mass for D = 2.5 and g = 0. 
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Figure 5: Dynamical fermion mass for D = 3.0 and g = 0. 
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Figure 6: Dynamical fermion mass for D = 3.5 and g = 0. 
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the case of the second order phase transition the critical temperature and 
chemical potential obey Eq. ([M|) . 

We introduce the eight-fermion coupling and evaluate the effective poten- 
tial. Since higher dimensional operators that we do not consider here have 
a larger contribution for a S> m-o, it is not valid to discuss such a parameter 
region with only the four- and eight-fermion interactions. We do not care 
about the stability of the effective potential and consider the parameter re- 
gion up to cr ~ 0(mo). Thus we can set the eight-fermion coupling g = —0.1 
or g = 0.1. 

The behavior of the dynamical fermion mass is shown as the chemical 
potential varies with a fixed temperature in Figs. [3 [8]and[9j A lower and a 
higher fermion mass is generated for a positive and a negative coupling, g, 
respectively. It is observed that the critical chemical potential is not modified 
by the eight-fermion coupling for the second order phase transition. As is 
shown in Fig.[7](b), there is only the first order phase transition for a negative 
coupling, g = —0.1, &t D = 2.5. In three dimensions we observe the first 
order phase transition at T = for 5 = 0. It becomes the second order one 
in Fig. El (a). 

In Figs. [7] (b), [8] (b) and [9] (b) the dynamical fermion mass at /i = is 
corresponds to the one predicted in Eq. (|21|) . Since it is a constant solution for 
the stationary condition of the effective potential, we observe a fiat behavior 
of the dynamical mass. The solution moves to the trivial one, md = 0, at 
the critical chemical potential through the first order phase transition for 
D = 2.5 . In Fig. [9] (b) we observe a small mass gap where the solution (j2ip 
drops to the other one ()19p then the broken symmetry is restored through 
the second order phase transition at the critical chemical potential. In three 
dimensions we observe both the behaviors. The constant solution moves to 
the trivial one at kT = and 0.2?7io. A small mass gap is observed for 
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Figure 7: Dynamical fermion mass for D = 2.5 as a function of the chemical 
potential /i with temperature kT/rriQ fixed at 0, 0.2, 0.4. 
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Figure 8: Dynamical fermion mass for D = 3.0 as a function of the chemical 
potential fi with temperature kT/rriQ fixed at 0, 0.2, 0.4. 
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Figure 9: Dynamical fermion mass for D = 3.5 as a function of the chemical 
potential /x with temperature kT/rtiQ fixed at 0, 0.2, 0.4. 



kT = 0.4mo. 

In Figs. [TOl [TT] and [T2] the behavior of the dynamical fermion mass is il- 
lustrated as a function of the temperature with the chemical potential fixed. 
The critical temperature coincides with the one in the four-fermion interac- 
tion model, 5 = 0, for the second order phase transition. We draw only two 
curves in Fig. [10] (a), since the dynamical fermion mass can not be generated 
at // = O.Smo for D = 2.5 and g = 0.1. It means that a positive g suppresses 
the chiral symmetry breaking and lowers the critical chemical potential for 
the first order phase transition. The dynamical fermion mass is given by 
Eq. dSU around T = in Figs. HO] (b) , E] (b) and 12] (b) . We observe a mass 
gap where the solution (fT8|) moves to the other ones, (fT7|) and (fT9]) . 

The transition between the solutions (|17p . (|18p and (|19p is clearly seen 
by observing the behavior of the effective potential. We show a typical one 
in Fig. [131 It should be noticed that the effective potential is normalized as 
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Figure 10: Dynamical fermion mass for D = 2.5 as a function of the tem- 
perature T with chemical potential n/rriQ fixed at 0, 0.4, 0.8. 
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Figure 11: Dynamical fermion mass for D = 3.0 as a function of the tem- 
perature T with chemical potential fi/rriQ fixed at 0, 0.4, 0.8. 



1.4 

1.2 

1 

I 0.8 

d* 0.6 

0.4 

0.2 









■ 




-^ 


%M=0 


1^= 


=0.4mo 


\\ 




li= 


0.8mo\ \\ 



1.4 

1.2 

1 

I 0.8 

^ 0.6 

0.4 

0.2 






0.2 0.4 0.6 0.8 1 1.2 

kT/iTio 



0.2 0.4 0.6 0.8 1 1.2 

kT/mo 



(a) 5 = 0.1. 



(b)5 



-0.1. 



Figure 12: Dynamical fermion mass for D = 3.5 as a function of the tem- 
perature T with chemical potential ^/niQ fixed at 0, 0.4, 0.8. 
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Figure 13: Behavior of the effective potential for g = —0.1, kT/rriQ 
0.2,0.3,0.4,0.5,0.6,0.7 and n/mo = 0.4. 



V{Q) = 0. For lower temperature the minimum of the effective potential is 
found at the solution (|18p . It moves to the trivial one, o" = 0, through the 
first order phase transition at the critical temperature in Fig. [13] (a). We 
observe that a local maximum appears between the solution (fTSj) and (fT9]l 
in Fig. [13] (b). Then the minimum of the effective potential continuously 
approaches to o" = 0. 



4.3 Phase boundary at finite temperature and ciiemical po- 
tential 

Evaluating the effective potential on the T — fi plane numerically, we find the 
phase boundary which divides the symmetric and the broken phase. We 
show the phase boundary for the four-fermion interaction model in Fig. [TJ] 
On the dashed lines the broken symmetry is restored through the second 
order phase transition, while the first order phase transition occurs on the 
solid lines. We observe the tricritical point where the first order critical line 
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Figure 14: Critical lines for g = 0. Figure 15: Critical lines for g = 0.1. 
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Figure 16: Critical lines for g = —0.1. 



meets the second order one for 2 < D < 3. 

In Figs. [15] and [16] we draw the phase boundary for g = 0.1 and —0.1, 
respectively. For the second order phase transition the critical temperature 
and chemical potential have no correction from the eight-fermion interaction 
in the leading order of the 1/A^ expansion. Thus we observe that the dashed 
lines are drawn on the same place, but some parts are overwritten by the 
solid lines. As is shown in Fig. [K] the same phase boundary is found for 
D > 3 when g is positive. In D = 2.5 the tricritical point is found at lower 
temperature. In Fig. [16] the first order critical line and the tricritical point 
appear even for D > 3. Only the first order phase transition is found at 
D = 2.5 for 5 = -0.1. 

5 Conclusion 

We have investigated the phase structure of the four- and eight-fermion in- 
teraction model at finite temperature and chemical potential in an arbitrary 
dimensions, 2 < D < 4. Evaluating the effective potential in the leading 
order of the 1/A^ expansion, we have calculated the dynamical fermion mass 
and the phase boundary. 

It is found that a lighter fermion mass is generated for a negative eight- 
fermion coupling, G2, (a positive g). Because of a constant solution (fT8]l for 
the stationary condition of the effective potential a new mass gap from the 
solution (HHD to ([TO]) appears for a positive G2, (a negative g). A similar 
behavior is found in Nambu-Jona-Lasinio model with 't Hooft and eight- 
quark interactions in four dimensions [19J. 

It is analytically shown that there is no correction to the second order 
critical line from the eight-fermion interaction. The first order critical line 
is changed and the tricritical point is modified by the eight-fermion interac- 
tion. Since the four-fermion interaction model is renormalizable below four- 
dimensions, only the irrelevant operator in the model is the eight-fermion 
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interaction. Thus the contribution from the irrelevant operator decreases 
compared with the renormahzable one near the second order critical line 
and disappears on the line. The same feature is found for the phase struc- 
ture in a spacetime with non-trivial topology [25j and a curved spacetime 
[231 [21]. 

The result seems to be inconsistent with the one in Refs. |16j and |26j . 
If we consider the model as a regularization of the one in four dimensions, 
the result depends on the regularization procedure. It is able to introduce 
a cut-off scale and define the model in four dimensions. Since the four- 
and eight-fermion interaction are irrelevant, both the operators have non- 
negligible contributes on the second order critical line. 

At the limit, /U — )• 0, the finite temperature phase transition is of the 
second order in the four-fermion interaction model. Thus the same critical 
temperature is obtained in the model with a negative G2. The first order 
phase transition can be realized for a positive G2. In this case a higher 
critical temperature is obtained. In the four-fermion interaction model the 
phase transition is of the first order only for 2 < D < 3 at the limit, T — )• 0. 
Thus the critical chemical potential takes the same value for a negative G2 
in 3 < D < 4. For a positive G2 the first order phase transition is found 
even in 3 < Z) < 4. 

In the present paper we only consider four- and eight-fermion interac- 
tions. The eight-fermion interaction produces a new local minimum for the 
effective potential and induces the first-order transition from the local mini- 
mum to the other. It is one of the characteristic features of the contribution 
from the eight-fermion interaction. Multi-fermion interactions which given 
by higher dimensional operators can also generate additional local mini- 
mums. Then multi-steps of the first order transitions induced in the phase 
diagram. 

We are interested in applying the results to some physical phenomena. 
The first order transition generates a gap for the energy-momentum tensor 
for fermion fields. It makes the deceleration parameter of our universe change 
suddenly. It also contributes the structure of dense stars. Thus we expect to 
observe such effect in cosmological and astrophysical objects at high temper- 
ature and/or high density. To discuss the physical phenomena we should fix 
the model parameters phenomenologically and consider some generalization 
of the model to be suitable for each energy scale. To investigate hadron 
physics we should include the flavor structure and current quark mass in the 
model. The attractive interaction between di-quark channel plays a decisive 
role to break the color symmetry at high density [3]. There is a possibility 
that super partners appear at SUSY scale. A super symmetric extension of 
the model is possible along the work in Ref. [27]. It may be important to 
study the dynamical symmetry breaking above the SUSY scale. We continue 
our work and hope to report on these problems. 
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